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Particles following circular trajectory widely exist in nature. Both lifeless circle particle and 
self-propelled circling microorganisms can be described by the same diffusion equation. We derive 
the diffusion equation of many circle swimmers from helical trajectory and Newton's law. The 
driven force includes both deterministic force and stochastic force. One special diffusion phenomena 
for Circle particle is that the density gradient and velocity gradient in one direction can induce a 
transverse flow in the corresponding perpendicular direction. This transverse flow will drive the 
particles to concentrate on a perpendicular direction to its original diffusion direction. Vortex in 
large scale is very easy to form due to the transverse flow. When many particles are trapped by the 
turbulent flow, there would appear a density bump upon the conventional exponential decay of non- 
circling particles. This non-equilibrium diffusion phenomena is originated from the mathematical 
circular trajectory itself, it preserves for different biological or physical system. An experimental 
evidence of this transverse flow showed up in a sperm diffusion experiment 
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I. INTRODUCTION 

In classical physics, there exist many familiar ex- 
amples for particles moving along closed orbital, such 
as planet in solar system, electrically charged particle 
in magnetic field. Many microorganisms swim along 
a helical trajectory until they reach a boundary, then 
they began circling near the boundary. The bacterium 
Escherichiacoli swims in clockwise circular motion near 
a soHd boundaryQQQidi- Many animal spermatozoa 
swims in circle when it is confined in a planar surface, the 
circling motion has a strongly preferred handedness Q ■ 
Artificial circle swimmers are constructed too, such as 
red blood cell attached by a chain of colloidal magnetic 
particles [§1, double faced nanorodes with different polar- 
ity on its two sides 0. 

Different self-propelled circle swimmers have different 
swimming mechanisms. The circling bacterium E.coli 
near a solid boundary is driven by a rotating helical flag- 
ella, theoretical study suggested that the hydrodynamic 
interaction between flagella and the boundary provide 
the driven force [2]. The sperm is also driven by flagella, 
but the beating pattern of the flagella is different from 
that of E.coli. The sperm flagella is theoretically mod- 
eled as left-right asymmetric beating ribbon, molecular 
motors inside the sperm trigger this beating [lo']. The 
over damped magnetic nanorod without flagella can also 
move in circles, however it is driven by a rotating ex- 
ternal field ^ll']. The dynamic mechanism for a circling 
rod particle swimming in confined channel is modeled 
as a Brownian circle swimmers [T2|. Both geometric and 
stochastic theor y a re proposed to study the motion of a 
single sperm[l3jfl4|. 

Recently an experiment was conducted to study the 
diffusion phenomena of millions of Arbacia punctulata 
spermatozoa in seawater, jelly coat solution and resact 
solution [l|. The experiment observation [l| is performed 
in a rectangular migration channel with sperm reservoir 
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and chemoattractant attached to its two ends(Fig. [l](a)). 
The diffusion image shows the sperm hken to concentrate 
one side of the channel(Fig. [l](b)). Sperms move faster 
near the boundary of the channeL 




(a) (b) 



FIG. 1: (a) The experiment device for the diffusion of sperms 
under attraction of chemoattractant in Ref. Ij. The length 
of the channel is 7000-9000 ^m, the width is 500 ^m, and 
the height of the channel is fOO fj,m. (b) A captured image 
of sperm diffusing in resact solution [J] . The bright dots are 
images of sperms. 

Motivated by this experiment, we developed a general 
theory on the diffusion phenomena of many circle parti- 
cles. We also derived the complete equation of motion for 
a single circle particle, one can apply this equation to de- 
scribe a Browian rod by including the length and the an- 
gle of rod. One can also apply this equation of motion to 
a single sperm by assuming an effective force of chemical 
stimulus. We shall start from the mathematical descrip- 
tion of a circling trajectory and Newton's law of motion. 
The effective driven force is of the most general formula- 
tion. One can set it as a deterministic force or a stochas- 
tic force for the convenience of physics, the formulation 
of equations are all the same. The diffusion equation for 
many circle swimmers is established on the conservation 
equation of total particle number. This diffusion equa- 
tion for the circles's center can be viewed as a member of 
the large family of Fokker-Plank equations, but the diffu- 
sion coefficients and drifting coefficients in our equation 
are not arbitrary functions, they are under strong con- 
strains and are determined by physical parameters. If 
we incorporate the local coordinates of the circle parti- 
cles into the diffusion equation, it becomes a sophisticate 
differential equation defined by two independent coordi- 
nation system, which is much more complicate than any 
conventional Fokker-Plank equations. We studied a sim- 
ple case of constant driven force. We found a transverse 
flow perpendicular to the density gradient. The direc- 
tion of this transverse flow depends on the direction of 
angular velocity. This transverse flow can help us to un- 
derstand the diffusion phenomena of Arbacia punctulata 
spermatozoa in seawater, jelly coat solution and resact 
solution [Tj- 



II. THE CENTRIPETAL FORCE EXPRESSED 
BY NONZERO ANGULAR VELOCITY 

This section is for the self-content of the manuscript, 
detail formulas are present to show that there is no as- 




FIG. 2: The coordination of one circling point-particle. 

sumption to introduce the centripedal force from a circu- 
lar trajectory, and no assumption is involved to separate 
the motion of the center and the motion of the particle 
itself. 

Newton's inertial law tell us, in the absence of force, a 
classical particle either is at rest or moves in a straight 
line with constant speed. If a particle draws a circu- 
lar trajectory, there must exist some effective centripetal 
force to drive it moving in that way. No matter how com- 
plex the mechanism might be to generated this force, the 
ultimate result of centripetal force is making a particle 
moving along circles. A nonzero angular velocity is a 
summary of the central pedal force. For specific physi- 
cal system, the expression of angular velocity depends on 
specific sources of interactions. In the following we will 
show how to express the centripetal force from a circular 
trajectory. 

We first consider a point particle moving along a circu- 
lar trajectory around the origin. The circular trajectory 
is described by r{t) = r^{t) + iry{t) = ro(t)e*"*. The 
coordination of the point-particle is 

x{t) — rQ{t) cos bJt, y{t) = rQ{t) sin cut, a;^ + = l.(f ) 

ro{t) is the radius of the circle. If ro{t) — aP,a > 0, 
we get an expanding spiral. If ro{t) ~ at^,a < 0, the 
trajectory is shrinking spiral. If ro{t) = ro + {asm{kt) + 
bcos{kt)), we get a fluctuating curve whose central line 
is a circle. We can also take ro{t) as a random variable, 
i.e., ro{t) = To -I- b£,{t)) with ^(t) as a random variable, 
then one will get a randomly fluctuating circle. For the 
simplest case, ro(t) = tq ~ const, the velocity of the 
particle is 

X — —tqIjJ svii{ujt) , y — +rQUj cos{u!t) , (2) 

one can verify the two components of the velocity, x and 
y satisfy 

x^+f^rluj^. (3) 

The two components of the velocity of a circling particle 
are not independent. The acceleration of the particle is 
given by the second derivative, 

X = —yuj, y = +xiu. (4) 

The motion of a circling particle along one direction pro- 
vides acceleration to its perpendicular direction as long 
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as w 7^ 0. We summarize Eq. ^ as, 



curve q(t) with respect to t give us the velocity, 



-mr X CO, CO 



(5) 



For the self-propel particle, such as bacteria or sperm, 
the centripetal force comes from interaction between the 
beating flagella and thick liquid where the bacteria lives. 
If the bacteria or sperm is in vacuum alone, no matter 
how hard its flagella is beating, it will not move along 
circles. The motion of a bacteria or sperm in vacuum 
alone depends on its initial velocity, if it is static in the 
beginning, it will not move at all even though its flag- 
ella is beating, if it has initial velocity, it will move in a 
straight line with constant speed, this is inertial motion 
which has nothing to do with the beating of flagella. For 
a circling electron, the effective driven potential could be 
magnetic field, angular frequency is given by w = eB/mc, 
B is magnetic field, e is electric charge. We put those de- 
tails aside, and summarize all those effective driven force 
into a non-zero angular velocity uj. 

If the center of the circle is not fixed and drifts around 
at low speed, the trajectory of the particle is a helical 
curve, it can be described by drifting circles. We in- 
troduce R{t) = Rx{t) + iRy{t) to express the instan- 
taneous center of the circle. The relative position of 
the particle with respect to the center is denoted as 
r(t) = Txit) + iry{t) = ro(t)e'"*, ro{t) is the relative 
distance between a particle and the center of circle. The 
complete position coordination of the circle swimmer is 



q{t) = x{t) + iy{t) ^ R{t) + r{t), 



(6) 



x{t) = Rx{t) + rx{t), r^it) =ro{t)cos{ujt), 
y{t) ^ Ry{t)+ry{t), ry{t) =ro{t)sm{Ljt). (7) 

This coordination provide us a convenient mathematical 
description for an arbitrary helical curve in two dimen- 
sional space. The shape of the curve depends on the func- 
tion of center coordination R{t) — Rx{t) +iRy{t) and the 
function of radius ro{t). For example, if the center accel- 
erates in X-direction with a constant acceleration a, and 
accelerates in Y-direction with acceleration 6, the center 
coordination functions are Rx{t) = ^at'^, Ry[t) = \bt'^. 
In the mean time, if we assume the radius of the circle 
is fluctuating periodically, ro(t) = -t- cos(fci)/100, the 
equation for a drifting helical curve is 



x{t) 



1 



-a r 



yit) = -be + 



ro 
ro + 



cos(fci) 



100 
cos(fci) 
100 



cos (jjt, 

sinojt. (8) 



From physics point of view, the tangent vector of the 
curve is the velocity of a particle which move along the 
curve. The acceleration can be expressed as the curvature 
of the curve. The time parameter t can be viewed as an 
arc length parameter of the curve. The first derivative of 



X — Rx{t) — vo{t) sin(wt) + fo(i) cos(ajt), 

y = Ry{t) + vo{t)cos{ujt) + ra{t)sm{ujt). (9) 

Rx{t) and Ry{t) are the drifting velocity of the cen- 
ter. vo{t) — ro(t)w is the speed of circling motion 
around the center. The mathematical definition of cir- 
cle swimmer requires that the circling speed of swimmer 
is much larger than the speed of its drifting circle cen- 
ter, i.e., vo{t) ^ Ra{t),a = x,y. For the opposite case, 
vo{t) ^ Ra{t), the trajectory of the particle is no longer 
drifting circles, and becomes an arbitrary curve. In that 
case, one can choose differential geometry to construct an 
equivalent theoretical description. The second derivative 
of q{t) with respect to t reads 

X = Rx{t) ~ yuj + ujRy{t) + ro{t) cos{ujt) 
— rQ(t)LO sm{Lut), 

y = Ry{t) + XLU — LdRx{t) + fQ{t) wiiuit) 

+ ro{t)ujcos{ujt). (10) 

The acceleration may be split into two parts: the accel- 
eration of the center, the acceleration of circling motion 
around the center. We multiply the mass on both sides 
of Eq. (fTO)) to demonstrate Newton's law. 



mx — m[Rx{t) + ro{t) cos{ujt)] 

— rri[yLU ~ uRylt) + rQ{t)uj sm{ujt)], 

my — m[i?j,(t) + fo(i) sin(ajt)] 

+ m[xuj — ujRx{t) + rQlt)!^ cos(iot)]. 

The absolute acceleration of the center is 

ax = [Rx{t) +ro{t)cos{ujt)], 
= [Ry{t) +ro{t)sm{ujt)], 
Vx = [x - Rx{t) +ro{t)cos{Lut)]. 



(11) 



[y - Ry{t) + ro{t) sin(wt)]. 



(12) 



The acceleration of a circling particle can be expressed 
into vector formulation, 



mz{t) = md — mV x uj, u 



(13) 



The non-zero angular frequency uj summarized the effec- 
tive centripetal force for the circling motion. The cen- 
tripetal force mV x cj is consumed by three motions, one 
is circular motion with respect to the center, the second 
consumer is the motion of the center, the third consumer 
is the expanding or shrinking of the radius, ma summa- 
rized the other external effective driven force Fextemai, 



external 



= ma, Fcentripetai = -mV X Q. (14) 



As long as we have a non-zero angular frequency uj, the 
effective driven forces from different sources are summa- 
rized into UJ. For an electron moving in magnetic field. 
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the non-zero w means the existence of magnetic field. 
For a self-propelled particle, such as bacteria or sperm, 
the non-zero lo summarized the effective centripetal force 
generated by the interaction between beating flagella and 
environment. 

III. THE DIFFUSION EQUATION OF 
CIRCLING PARTICLES IN TWO DIMENSIONS 

A. random fluctuation upon deterministic 
trajectory 

Circling motion is a kind of deterministic dynamics 
driven by a deterministic force. When a particle is 
surrounded by many thermally fluctuating molecules, 
the constant collision between the particle and other 
molecules would cause stochastic deviations from the de- 
terministic path. This stochastic fluctuations is driven 
by a random noise force. The amplitude of stochas- 
tic fluctuation £,x{t) should be much smaller than the 
radii of the circle so that people can observe a fluctu- 
ating circle trajectory. If the fluctuating amplitude is 
much larger than the radii of the circle, one would just 
see some random trajectory instead of circular trajectory. 
So we take the deterministic path as dominant trajectory, 
and introduce a small amplitude of random fluctuation 

x{t) = Rx{t)+ro{t)cos{Lot)+il{t), 

y{t) = Ry{t)+ra{t)sixi{ujt)+ey{t). (15) 

The random fluctuation f (i) can be viewed as the fluc- 
tuation of the center coordination, 

x{t) = [Rx{t)+il(t)]+r^{t)cos{ujt), 
y{t) = [Ry{t)+ey{t)]+ro{t)^H'^t), 

(16) 

^{t) can also be taken as the fluctuation of the radii of 
the circle, 

x{t) = Rx{t) + [rQ{t)+(i{t)]cos{ujt), 
y{t) = Ry{t) + [ro{t)+ey{t)]sin{u;t). 

^x{t) -^roit), ^y{t)^ra{t). (17) 

The angular frequency is not always constant for some 
self-propel particles. So we can introduce a random vari- 
able ^(t) into the deterministic angular frequency, 

x{t) = R4t)+roit)cos{[ujit)+^4t)]t), 
y{t) = Ry{t)+ro{t)s\ii{[u:{t)+^y{t)]t). 

(18) 

One can choose different coordination for different physi- 
cal systems. If one can fix the center and the radii of the 
circle at the same time, the fluctuation may be combined 
into the angular frequency. If the center and radius of the 



circle is not fixed, for example, a randomly swimming 
bacteria, it is more convenient to summarize the small 
stochastic motion by fluctuating center or fluctuating ra- 
dius. For the most complex case, we can combine these 
three types of fluctuation together. The random fluctu- 
ation of the center coordinates and angular velocity are 
independent. One can choose a convenient description 
for different cases. 

We shall consider the diffusion phenomena of millions 
of circles in the following, each of them obeys an equation 
of motion, 

X = Rx{t) — wo(i) sin(wt) + fo{t) cos{ujt), 

y = Ry{t) + vo{t)cos{ujt) + ro{t)sin{ujt). (19) 

where uj — uj + (t) , the angular frequency is fluctuating 
around a central frequency. The solution of Eq. ([TO)) is 
a drifting circles. 



B. Diffusion equation of circle swimmers in liquid 
with low Reynold number 

When we talk about diffusion phenomena, it usually 
means the evolution of density distribution for a large 
amount of particles in an area with or without boundary. 
Density is defined by the number of particles in one unit 
area. Circle particle has two microscopic scales, one is 
the size of the particle itself, the other is the size of the 
circle. We can define density as the number of particles 
in an unit scale of the particle itself. If we fix the center of 
the circle, the density function for one particle is a sharp 
peak circling along circular trajectory. The average of 
these peaks in a longer time scale just sit at the center 
of the circle. We will not see diffusion phenomena unless 
the center of the circle moves. So we should choose an 
appropriate unit scale to get an appropriate definition of 
density function for diffusion. The unit spatial scale is the 
radius of the circle, the unit time scale is the periodicity 
of the circling motion, 

Ax«2ro, Ay^2ro, Ai w Tq, (20) 

where Tq measure how much time is spent for a particle to 
draw a circle. We can either project the mass of particle 
onto the center of the circle, or fill it homogeneously in 
the unit disk 7rr§. We assume that several particles are 
allowed to swim along the same circle. This is to meet the 
experiment observation of sperm swimming that two or 
more sperms can swim together 6] . We assume there can 
be at most s particles in one circle. The biggest number 
of circles in a finite area 5 is / = S/tttq. We can define 
density as the number of particles in one circle. 

If the step size of experiment measurement is much 
larger than the radius of the circle, we can choose a larger 
unit scale 

Ax > 2ro, Ay > 2ro, At > Tq, (21) 
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The density is defined by the number of particles in unit 
box AxAyAt. In an thermal environment, the phases 
of circling motion for different particles are usually ran- 
domly distributed. The average of these relative circling 
velocity with respect to the centers of the circles is zero. 
The internal motion of the particle inside the unit box 
can be neglect for diffusion phenomena. However the 
centripetal force induced by the circular motion will be 
consumed by the center. 

We separate the motion between the circling and the 
drifting of centers, the circling motion is described by a 
deterministic equation 



'VQ{t) sm{ujt) + fo(i) cos(wt), 
-I'oit) cos(w<) + ro{t) sin(a;i). 



(22) 



The drifting motion of the center is determined by dif- 
fusion dynamics which does not concern the motion of 
particle within the scales of a circle. In the diffusive dy- 
namics, we only need to study the dynamics of the center. 



Rx{t), Ry{t)- 



(23) 



We can combine the motions of the center and the small 
circling motion together to get the equation of motion 
for every particle, but for diffusion phenomena, it will 
not make much difference. 

When a swimmer is placed into a collection of thou- 
sands of swimmers which has inhomogeneous density dis- 
tribution, it has larger probability of entering into the 
low density region, but less probability of moving into 
the high density region. Because there are more free 
spaces in low density region, and less free space in the 
high density region. The density gradient acts as an ef- 
fective driven force to drive a swimmer from high density 
region to low density region. We denote it as diffusive 
driven force Fdif fusion oc -{SIrNs)/Ns. 

Physical circle swimmers have a finite body volume 
and an effective circling volume. If a swimmer tries to 
push itself into a crowd of swimmers, it has to expel 
the molecules and other swimmers in that area, it will 
inevitably encounter friction resistance. The faster it 
moves forward, the stronger resistance it will receive, so 
we take the resistance as a negative force proportional to 
velocity , Ff„cUon oc -r]V. 

According to equation of motion ([T^ and p^ . we add 
the external driven force F^xternai and centripetal force 
Fcentripetai upou the Center. We are dealing with millions 
of particles, one reasonable choice of elementary unit is 
to take a small unit box of many particles. We neglect 
the internal motion of individual particle inside the box. 
The unit box is viewed as a whole to interact with other 
boxes, the equation of motion for a square of particles is 
given by Newton's law, 

dV ^ ^ 

mNs — = -mNsV xu + N.F-DVNs-mNsT^V. (24) 
ot 

•q indicates the friction coefficient. In fact, if we divide 
both the two sides of Eq. with the density iVg, it 



leads to an effective equation of motion for single particle. 



mR = ~mR x co + F — D- 



mriR. (25) 



We can equally split the unit box's centripetal force, ex- 
ternal force F and friction force mNsTjV into TVs pieces, 
each particle inside share one piece, then we will get Eq. 
((25)) for single particle in the absence of a special term, 
the effective diffusion force 



fdiffus 



-D 



(26) 



Density gradient is a conception of many particles, for 
a single particle, density gradient does not exist. This 
diffusion force term only exist when the particle is placed 
in a many particle system. The external driven force F 
may include all different kinds of forces, deterministic or 
stochastic, i.e.. 



Fx = f^-dRxV(R) 

Fy = f^-dRyV{R) 



m 

m 



(27) 



V{R) is a trapping potential, is a deterministic driven 
force, /^(t) = rn£^x{t) and fy{t) = m^y{t) is white noise 

random force. If we identify F as white noise random 
force, Eq. (P5|) becomes a generalized Langevin equation 
with a special diffusion driven force f diffusion- Equation 
of motion ([25)1 includes two coupled equations. 



mRx = —mujRy — D. 



dR.xNs 



mriRx + Fx, 



mky = muR,-Dy^^-m^Ry + Fy, 

^ s 

Here the diffusion coefficient Dx and Dy are state depen- 
dent function, Dx = Dx{x,y) and Dy — Dy(x,y) 

We assume the particle swims in a thick liquid with 
low Reynold number, the friction on the center is much 
larger than the inertial term of the center, the motion of 
the center is over damped. Then we set the inertial term 
as zero. 



mRx = 0, mRy 



0. 



(28) 



Noticing here it is the center's inertial term, the motion 
of the particle itself is not over damped since it has self- 
driven force. This self-driven force comes from the in- 
trinsic beating pattern of the flagella. They were born to 
swim in circles. In the smaller scale inside the circle, the 
particle generates an effective driven force by its internal 
motors, this effective driven force canceled the friction 
force along the circle, in the meantime generated an ef- 
fective centripetal force to make itself follow a circular 
trajectory. Part of this effective driven force is consumed 
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by the center through the term —mR x ui. The formula- 
tion of velocity is derived by solving Eq. 



Rx 
Rv 



-A, 



D 



vv 



dB,yNs 



^yyFy ^yx^x- 

(29) 



where the diffusion tensor D and drifting tensor a are 



DxV ^ _ DyUJ^ 



D - ^ D 
V 



'xy 



mA^ ' 



A^' 



(30) 



A^ = rj^ + bj^ is a normalization factor. We express Eq. 



29)1 by a brief formula, 



i?,; — 



aijFj{x,y) 



3 



(31) 



Combining the over damped center equation and the de- 
terministic circular motion, we get the complete equation 
of motion 

^ — ^xy _ ^xx ~^ ^xx^x ^xy^yt 

— ro(t)LU sin(ti;<) -|- ra(t) cos(a;i). 



dRxNs 



dRyNs 



^yyFy '^yxFx 



^yX J^yy ^ -T Uyy± y ^ Uyx-L X- 

+ ro{t)u! cos{ujt) + fo{t) sm{ujt). (32) 

The solution of this equation is a drifting circle. This is 
the dynamics of a single particle. When we considering 
the diffusion phenomena of millions of particles, we care 
about how many particles is distributed in certain area 
at time t. The unit scale used to measure the diffusion is 
usually larger than the radius of the circle. The circling 
phase of the particles inside the unit box are random, so 
we can neglect the deterministic circling part, and only 
take into account of the motion of the center. The total 
number of swimmers flow into an unit box must equal 
to the total number of particles flowing out of the unit 
box. The particle flow is Jq = NgV. The differential 
formulation of the mass conservation equation of particles 
is 



dtN, +Vr-Jo^O. 



(33) 



Substituting Eq. (1^ into the mass conservation equa- 
tion, we derived the diffusion equation. 



dt 



E 



dRt [a-tjFj (x, y)Ns - Dij {x, y)dRjNs 



0. 



(34) 



The friction coefficient ry may be state dependent when 
the particle is in a complex environment, so we introduce 
the friction function rj = ri{x,y) to meet more general 
cases. The drift tensor a = (^{x, y) is state dependent as 
well. We introduce the drift function. 



Ti{x,y) = ^aij{x,y)Fj{x,y), 



(35) 



to express the diffusion equation p4p into the formulation 
of Fokker-Planck equation, 

^ + J2[dmmx,y)N,] - dR,[D,j{x,y)dR,Ns]] = 0. 

ij 

(36) 

If the diffusion tensor is homogeneous in space, they 
are constant without any dependence of space time vari- 
ables. Then we get a simpler diffusion equation, its ex- 
plicit formulation is. 



dNs 
dt 



~ [Dxy + Dyx]dRxdRyNs 

- DxxdlxNs-DyydlyNs 

+ [^yyPy + ^yxFx] dRyNg 

+ [o'xxFx + f^xyFy] dR^Ng 

+ [^yy^RyFy -f (JyxdRyFx] 

+ [c^xxdRxFx + OxydRxFy] Ns 

= 0. (37) 



One simplification of this complex diffusion equation is 
to set the diffusion coefficient along X-direction and Y- 
direction as the equal number, Dx — Dy = D. Then the 
diffusion tensor becomes simpler. 



Dr = Dx 



D 



yy 



D-q Dlu 

r, Dyx - -Dxy - (38) 



7nA^ ' 



C. Transverse force flow induced by velocity 
gradient 

One special character of circle particle is that the two 
components of its velocity are not independent. If its av- 
erage velocity increase in X-direction, its average speed in 
Y-direction will increase as well so that it can draw a cir- 
cle, otherwise, its trajectory just becomes ellipse, or any 
other irregular geometry. Thus if there exist velocity gra- 
dient, a transverse flow will be induced on the boundary 
between fast swimmers and slow swimmers. The veloc- 
ity field for circle particle is expressed by density function 
and external driven field. We define the transverse flow 
induced by velocity gradient as 



J 



(39) 



where uj = we^ for two dimensions, denotes the nor- 
mal vector oi X — Y plane, w < means the swimmer 
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is rotating in clockwise direction, w > means the cir- 
cle swimmer rotates in counterclockwise direction. In 
physics, the derivative of velocity is acceleration. So the 
vector flow J is actually the effective driven force field. 
For the velocity field Eq. (|29|) . the transverse force flow 
is 

Jy — UjdnxWxxFx] + ^diixWxyFy], 

+ ujdRx[Da;ydRylnNs] + ujdRx[DxxdB.x'^T^Ns] 

Jx = UjdB.y[<7yyFy] + UjdRy[ayxFx]. 

+ UjdRy[DyxdB.X^riNs] + LL!dR,y[DyydRy\nNs].{4Q) 

When the circle swimmer switch its circling handedness, 
i.e., from clockwise to counterclockwise, or vice verse, the 
transverse flow switch its direction as well. The switch- 
ing of circling handedness is simply flipping sign of the 
angular velocity, uj — > — or —uj — )■ uj. According to 
Eq. pp]) . the direction of flow Jy will flip to — Jy, so does 
Jx ■ If the velocity gradient only exist in X-direction, the 
transverse flow J is along Y-axis(Fig. [3] (a)), —uj and uj 
leads the transverse force flow to opposite direction. We 
can predict an observable experiment phenomena accord- 
ing to the flow fleld equation J = (Jx, Jy). We mix two 
groups of circle swimmers with opposite circling handed- 
ness together, and let them diffuse in a rectangular chan- 
nel. The particles will splits into two separated groups 
during diffusion. All the particles in one group have the 
same circling handedness. The two groups are moving in 
opposite direction, one goes to positive Y-direction, the 
other goes to negative Y-direction (Fig 13] (b)). 



of the same size as one circle. We can classify this small 
vortex by vorticity which is defined by 



Si 




Ai. 






o 


o 

►J 




o 





op oo 



(a) 



(b) 



FIG. 3: (a) The transverse flow induced by velocity gradi- 
ent, (b) Swimmers with opposite circling handedness will 
split during diffusion. 



D. Vortex generated during diffusion 

The sea urchins spermatozoa confined in planar sur- 
face are randomly swimming at low density. When the 
sperm density is increased to about 2500 spermatozoa 
per mrn? , spermatozoa self-organized into dynamic vor- 
tex array with local hexagonal order Q. Every vortex 
contains 10 ± 2 sperms swimming with the same hand- 
edness. This kind of vortex is the smallest vortex. It is 



where the velocity field is 

rx — —rQUJsiii{ujt) = —ujry, 
fy = +r{)Uj cos{ujt) ~ ujrx- 



2uj. 



(41) 



(42) 



This vorticity F^ quantifies the vortex observed in sperm 
experiment's] . If we go to larger scale, the non-zero an- 
gular velocity upon the centers of those small circles will 
also generate vortex in larger scale. The vorticity of the 
center's velocity field is. 



Fv 



\/xV = dRxRy - dRyR 



(43) 



here the velocity field of the center is given by Eq. (|29|) . 
We consider the special case that the diffusion tensor 
are state- independent, and the diffusion coefficient in X- 
direction equals the diffusion coefficient in Y-direction, 
Dx — Dy, then the diffusion tensor and drifting tensor 
satisfy the relations. 



J^xx — F^VV: F)xy — Dyx, 

CTxX = CFyy, (Jxy — —<Jyx- (44) 

Substituting Eq. into Eq. gS]), we get 

= CFxx(dRxFy — dRyFx) + CFyx{dRxFx + djlyFy) 

+ Dxy{dRydRy\nNs ^ dRxdRx\nNs). (45) 
According to the Green-function theory, 

dRydRylnNs + dRxdRxlnN, = 2^5{Ns), (46) 
we can simplify Ty further, 

Ey = CTxxidRxFy — dRyFx) + Cryx{dRxFx + dRyFy) 

+ Dxy27rS{Ns), (47) 

The particles are circling around a common center, they 
are pulled away from the center by centrifugal force, so 
the density vanishes at the center. We can technically 
split the density function into two components, and ex- 
press the density function as Ns — \/ N"^ + iV^ . Then we 
can identify every vortex with a topological number by 
applying Duan's topological current theory 0. It should 
be noticed that there are many vortex field here. 

More over, the transverse force flow induced by velocity 
gradient may also have vortex configuration. We measure 
this kind of vortex by 



VX J — dRxJy — dRyJx 



(48) 



Fj defines a higher order of vortex. Both Fj and Fy are 
determined by the density distribution and external field. 
When the density distribution evolute following diffusion 
equation (j62|l . the vortex configuration will change cor- 
respondingly. The vortex generated by the external field 
is controlled by external potential. 
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E. The fluctuation of diff"usion coefficient 

The diffusion tensor depends on friction coefficient, dif- 
fusion coefficient and angular velocity. The angular ve- 
locity of a circling particle confined in two dimensions 
has only one component uj^- Bacteria swim in circles 
near solid surface. Experiment found the angular fre- 
quency of bacteria's circling has strong fluctuations [l^. 
The angular frequency actually measure the strength of 
centripetal force. If we fix the center of the circle, A fluc- 
tuating angular frequency just make the particle move 
faster or slower along the circle, the center is static. If 
we release the center, it will consumed part of the cen- 
tripetal force. So the diffusion coefficient obtained above 
strongly depend on the fluctuations. We introduce a fluc- 
tuating vector ^z{t) into the angular velocity vector 



(49) 



For the state independent diffusion coefficient, the diffu- 
sion tensor with angular fluctuation is. 



DxV 



(50) 



If the fluctuation reduced the angular velocity, both dif- 
fusion coefficient D^x and D^y will increase(Fig. 14]), but 
the Dyx will reach climax first(Fig. [5]). If the fluctua- 
tion increase the angular velocity, both the diffusion co- 
efficients Dxx and Dxy will decrease, but Dxy decreases 
much slower than Dxx ■ If we take the ratio between Dxx 
and Dt 



(51) 



one can see clearly that the increase of frequency will 
make the off-diagonal diffusion stronger. As the total 
number of particle must be conserved, the increase of 
off-diagonal diffusion will decrease the diagonal diffusion. 
When the friction coefficient 77 is state dependent, i.e., 
J] = r]{x,y,t), the behavior of diffusion coefficient is de- 
termined by the complete coefficient Eq. (|50l) . 



IV. THE DIFFUSION EQUATION FOR = Dy 

We first study a special case of diffusion equation 
([M]). The diffusion tensor and drifting tensor are state- 
independent. We take the diffusion coefficient Dx as an 
equal value of Dy , 



Dx=Dy= D. 

Then the off-diagonal diffusion coefficients and drifting 
coefficients are. 



-Dx 



Duj 



(53) 




FIG. 4: The diagonal diffusion coefficients (Dxx^Dxy) under 



random fluctuation ^. We have set rj/m 
frequency is cj = 2. 



1, and central 




FIG. 5: The off-diagonal diffusion coefficients {Dxy~—Dyx) 
as a function of random fluctuation 5. Parameters are given 
by rj/m — 1, lj — 2. 

The external driven force is constant both in X-direction 
and in Y-direction, 



^x — Cx 



(54) 



If ^ 0, the particle would have a drifting velocity in 
Y-direction. For a non-circling particle, if Fy = 0, the 
density center will not drift in Y-direction. But for circle 
particles, even if Fy = 0, the density center still has a 
drifting velocity in Y-direction. To see this more clearly, 
we take a further step of simplification, 



Fx^Cx^Q, Fy^O. 



(55) 



There is no external force in Y-direction now. For the 
circle particle, the angular frequency uj is not zero. If 
the particle is circling is in counterclockwise direction, 
a; > 0. If the circling is in clockwise direction, w < 0. 
The diffusion equation ([M)) now becomes. 



dt 



^D, 



,dl,N, + D 



yy^Ry^s 



b diiyNs - a OrxNs 



(56) 

(52) The drifting velocity is a = ctxxFx, b ~ cTyxFx- The 
drift of density center along Y-direction is driven by off- 
diagonal diffusion tensor ayx- Since there is no external 
driven force in Y-direction, this drifting velocity is a spe- 
cial character of circle particles. 

The analytic solution of diffusion Eg. (1561) with infinite 
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boundary condition is 
1 



exp[ 



(i?^ - atf {Ry - htf 



ADrt 



ADrt 



(57) 



The center of sperm density is drifting under the action 
of external driven force. We integrated the density along 
Y-axis and derived the density evolution along X-axis, 
Ns{x, t) — Ngdy, where h is the width of channel along 
Y-direction. If the drifting velocity is much larger than 
the diffusion coefficient, we would observe apparent drift- 
ing of the density center(as shown in Fig. [S] (a)). When 
the diffusion coefficient is much larger than the drifting, 
one would see a rapid decay in beginning, then there 
appears a plateau due to the drifting(Fig. H] (b)). We 
integrated X out to keep the Ng as a function of {y,t), 

Ns{y,t) = jj* Nsdx. One can see the drifting velocity 
of density center in Y-direction is b. When particles are 
diffusing in X-direction, their density distribution in Y- 
direction is not homogeneous, they concentrates to either 
positive Y-direction or negative Y-direction. 




FIG. 6: (a) The density evolution of solution equation (|57p 
for a/Dr — 7. The drifting velocity is larger than diffusion 
velocity, so one can observe a moving pump, (b) The density 
evolution of solution equation (|57p for a/Dr = 1/4. The 
drifting speed is slower than the diffusion speed. 

There exist a typical vortex configuration for the sim- 
ple case of = Dyy = D,. (Fig. [7]). For the analytic 
solution of diffusion equation (1571) , we can get an explicit 
formulation of the velocity field Vr — (14, Vy). The vor- 
ticity of this vortex is given by the vorticity Eq. (|43|. 
Substituting the density function Eq. (1571) into the ve- 
locity Eq. (|29p . then we can express the vorticity Eq. 
(|43| explicitly as 



Tv = D^y[dRydRyL{R) + 8 R^O R^R)] . (58) 

where the function L{R) is 



L{R) = 



(i?, - atf {Ry - btf 



ADrt 



ADrt 



-A-KDrt 



(59) 



The vorticity of the transverse force flow J also has ex- 
plicit expression 



Tj = D,yUj[dl,dRyL{R) ~ dlydR,L{R)] 



(60) 



as 



'III 



FIG. 7: The velocity field V = (14, K/) for density solution 
(|57p at time t = 1. The drifting parameters are a=l, b=l. 
Diffusion coefficient Dr = 1. 



where the function L{R) is given by Eq. (|59p . Fj quan- 
tifies the vortex in the space of higher order derivative of 
density distribution. 

If the particles are confined in a finite rectangular re- 
gion, the solution of diffusion Eq. ((5B)) do not have ex- 
plicit expression. We have to do numerical computa- 
tion to study the density evolution Ns{x,t) = Nsdy. 
We take a rectangular region covered by < x < 10, 
^ y ^ 5. The diagonal diffusion coefficient are 
Dxx — Dyy ~ Dr = 1. The density evolution for a 
small drifting velocity a = & = 1 is an exponential decay 
in the beginning (Fig. |H|). As the diffusion goes on, the 
exponential decaying curve is approaching to a straight 
line. Further more, we computed the density evolution 
for larger drift velocity a = b = 18. The density distribu- 
tion curve first transform from an exponential decaying 
curve to a curve with plateau, and then approaches to a 
linear distribution (Fig. 




FIG. 8: The numerical evolution Ns{x,t) = Nsdy of dif- 
fusion equation ((56} in a rectangular region. The drifting 
velocity is small a=b=l. 
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FIG. 9: The density evolution Ns{x,t) — Nsdy of diffusion 
equation ((56} in a rectangular region. The drifting speed is 
very large, a=b=18. 

V. THE DIFFUSION EQUATION FOR / Dy 

When Dx ^ Dy, the diffusion equation will include a 
cross-diffusion term, DyxdnxdnyNs. If the angular veloc- 
ity of the particle is zero, = 0, this cross-diffusion term 
vanishes. The density distribution for uj > and a; < 
are different. The numerical evolution of diffusion equa- 
tion shows, the density moves to opposite directions with 
respect to w > and uj < 0. We study a special case as 
an example to demonstrate the difference between uj > 
and w < 0. The diffusion coefficients and drifting tensor 
are state independent. The external force is constant, 

Fx — Cx, Fy = Cy, (61) 

The diffusion equation for D^^Dy is, 

^ = DxxdlxNs + DyydlyN, 

+ {Dyx + Dxy)dRxdRyNs 

- b dRyN, - a DrxNs. (62) 
The drifting velocity is (J — ^xxCx ^xyCy, h — ^yy^y ~^ 

^yx Cx • 

A. The diffusion equation for uj > Q 

We first study the diffusion phenomena of a positive 
angular frequency, w > 0. Diffusion Eq. ([62)) has no an- 
alytic solution. We performed a numerical computation 
of the density evolution in the two dimensional plane of 
Rx — Ry. In the beginning, the particles are distributed 
along Y-axes at i?a; = 0, the initial density is Ng — 10. 
The diffusion coefficients are Dxx = 0.9 and Dyy — 1. 
We choose a large drifting velocity, 6 = 5 and a — 15. 
The off-diagonal diffusion coefficient is Dyx + Dxy = — 1 
which corresponds to a positive angular frequency. One 
snapshot of the density evolution at i = 15 is shown in 
Fig. [TO] (a), the light color demonstrates the high den- 
sity, the dark color shows the low density. The center 



of high density region shifts to positive i?j/-direction(Fig. 
[ini(a)). 

Later we increase the off-diagonal diffusion tensor and 
drift tensor to Dyx + Dxy — —2, b=10, but the diagonal 
diffusion tensor remains invariant. A snapshot of density 
evolution at t = 15 shows the drift in Y-dircction speeds 
up(Fig. [TO](b)). In X-direction, the density center does 
not reach as far as the diffusion in low angular frequency. 




(a) (b) 

FIG. 10: (a) The density distribution Ns{Rx, Ry) at time 
t — 15 for the numerical evolution of diffusion equation (I62|) . 
The bright color indicates high density, dark color indicates 
low density. The particles concentrate on the line Rx = at 
time f = 0. parameters are Dyx + Dxy = —1, oi > 0, fo = 5 
and a = 15 . (b) The density distribution Ns{Rx, Ry) at time 
t = 15 for Dyx + Dxy = -2, 6 = 10, oj > and a = 15. 

The vortex field for Dx ^ Dy is very complicate. We 
computed the vortex field of diffusion Eq. (|62p with re- 
spect to the density distribution Fig. [10] (a). All the 
parameters have the same value as that of Fig. [TO] (a) . A 
vortex line(Fig. [TT] (a)) exist along the high density re- 
gion of Fig. [TO] (a). Local velocity vectors fiow out of the 
vortex line. The vortex line is a source of particles. This 
is just a vectorized demonstration of diffusion. Diffusion 
is a phenomena of many particles flowing from high den- 
sity region to free space. As time goes on, the center of 
the vortex is moving along with the density center dur- 
ing diffusion. At time goes on, the center of the vortex is 
drifting with a velocity determined by drifting parame- 
ter a and b. The effective radius of the vortex line is also 
expanding in diffusion. We computed the corresponding 
vorticity function(Fig. [TT](b)). The high vorticity region 
is covered by bright color. The bright islands of high vor- 
ticity are situated at the center of vortex lines(Fig. [TT] 
(b)). High vorticity islands are distributed in the upper 
half plane along Y-axes. 

The vectors in Fig. [TT] (a) is the velocity field of parti- 
cles. The direction of the vector is exactly the direction 
of particle's velocity at time i = 15 at certain point in 
Rx — Ry plane. The vortex line is a source of particles. 
Particle are fiowing out of the vortex region. Particles 
are trapped in the turbulence region, this will inevitably 
increase the density at that region. In this sense, the 
vorticity distribution is also the distribution of density 
center. When particles are expanding from high density 
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(a) 




FIG. 11: (a) The velocity field V = (14,14) corresponds to 
the density distribution of Fig. [10] (a), w > 0. (b) The cor- 
responding vorticity distribution for the vortex configuration 
of the left panel. The light color indicates high vorticity, dark 
color indicates low vorticity. This figure is a snapshot at time 
/; = 15 for the positive angular frequency case, cj > 0. 

to free space, the radius of vortex is also expanding. The 
region covered by non zero vorticity becomes larger as 
diffusion goes on. 

The evolution of the velocity at next step is determined 
by the acceleration of the particles. For circle particle, 
we can see where the particles is going to move from the 
transverse force flow J — {J^, Jy). The transverse force 
flow J is determined by density distribution, 

Jy = DxyUjdRxdRyXnNs + Dx^ujdRxdRj\nNs, 

Jx = DyxUjdRydRx'XnNs + DyyUjdRydRylnNs. (63) 

The transverse flow distribution is constantly changing 
from one configuration to another during diffusion pro- 
cess. We computed the corresponding transverse flow 
field J — {Jx,Jy) (Fig. [12]) with respect to the density 
distribution of Fig. [TU] (a). The flow field indicates the 
acceleration direction of particles. Fig. [T^ shows the par- 
ticle on the upper boundary flows straightly to positive 



X-direction, the particles in the lower half plane flows up 
to the upper half plane. In the high density region, there 
exist a turbulence region, vortex configuration come into 
being there. Compare the flow field (Fig. [T^ with the 
corresponding density function (Fig. [10] (a)), we see that 
the turbulence region situated right at the high density 
region. 




FIG. 12: The effective force flow field J = {Jx,Jy) cor: 
sponding to the density distribution of Fig. [10] (a), oj > 0. 



B. The diffusion phenomena for cj < under the 
same parameter setting as cj > 

In last section. We have studied the diffusion equation 
for oj > 0. The angular velocity in the diffusion Eq. ([5^ 
is positive, that means the circle particles are circling in 
counterclockwise direction. If the angular velocity switch 
from +0J to —OJ, the off-diagonal diffusion tensor and drift 
tensor will switch their sign simultaneously, i.e.. 



Jyx 



'yx- 



(64) 



So the off-diagonal diffusion coefficient of cross diffusion 
term will flip a sign, 



{D. 



yx 



D 



xyj 



(^Dyx ~t~ Dxy\ 



(65) 



The drift velocity of the density center will switch to 
negative Y-direction. The diffusion in X-direction is still 
in positive X-direction. For {Dyx + D^y) = 1, b=-5, 



a=15, we performed the numerical evolution of Eq. ()62|) . 
the density distribution at time f = 15 is shown in Fig. 
[T51 The particle concentrated on lower half plane of Y- 
direction. 

We computed the corresponding velocity field(Fig. [14] 
(a)) with respect to the density distribution of Fig. [TS] 
one can see the vortex line shifted the lower half plane 
too, so does the turbulence region(Fig. [14] (a)). Numer- 
ical evolution of the corresponding vorticity distribution 
also shifted to the lower half-plane (Fig. [14] (b)). The 
corresponding transverse force flow J = {Jx,Jy) for the 
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FIG. 13: One snapshot of the numerical evolution of density 
distribution of equation (|62l) at time t = 15, < 0. The 
bright color indicates high density, dark color indicates low 
density. 

density distribution Fig. [13] shows the particle flow in 
the upper half plane now point down to the lower half 
plane(Fig. 1151) . The most turbulent region now shifted 
to the lower half plane. The vortex sits right at the high- 
est density region. Now we can come to a conclusion 
that, when many circle particles are diffusing across a 
rectangular channel along X-axis, circles with positive 
angular frequency drift to upper boundary(positive Y- 
direction), while circles with negative angular frequency 
drift to bottom boundary (negative Y-direction). 
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VI. DIFFUSION IN A CYLINDRICAL DENSITY 
GRADIENT 

If initial density distribution of particles has cylindrical 
symmetry, it is more convenient to express the diffusion 
equation ([Ml) by polar coordination system. Suppose the 
particles concentrate in a small disk with the origin as its 
center, then the gradient of density would has cylindrical 
symmetry. Newton's equation for circle's center in polar 
coordination is of the same form as that in Cartesian 
coordination, 

m-^ = —mV X uj + F — D ^ — mrjV , (66) 

Here the velocity \sV — (Vr, Vg). The two components 
of the velocity are 

Vr^R, Vg^RB, VR^Vg. (67) 

Newton's law give us the same equation of motion, 

mVe = -muVR 77—77 m'qVg + tg, 

R J-^s 

uiVr = mulVe - -D_r-|^ - mr]VR + Fr, 



FIG. 14: (a) The velocity field V = iVa:,Vy) for uj < 0. (b) 
The corresponding vorticity distribution Ty of diffusion equa- 
tion ([621) for UJ <0. 




2 4 6 8 

FIG. 15: The transverse force flow fleld J = {Jx, Jy) for uj < 
with respect to the density distribution of Fig. 1131 All the 
parameters are the same as that for Fig. 1131 
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(a) 



(b) 



FIG. 16: (a) The difFusion of circling particles across a cylin- 
drical density gradient, (b) The diffusion phenomena with 
two blocking walls situated at = and 6 = n. 



Dg is the difFusion coefficient along the circle. Dji is the 
diffusion coefficient along the radius. In the over damped 
case, we set mVe = and mVn = 0. Then we can solve 
the equation above to get explicit equation of velocity, 



Re ^ -D, 



OR- 



R 



-Due 



1 dgN, 

R Ns 



Dg 



1 dgNs 

R 

dRN, 



^Fg + agRFu, 



Here the difFusion tensor D and drifting tensor a are 



Dgg = 

Drr = 

<^Be = 

ctrr. = 



Drt] 



DgR = - 
Drb = 



DgUJ 



croR = 



mA^ ' 

UJ 



(69) 



rj^ + is a normalization factor. The equa- 



where A^ 

tion of a circle in polar coordination is 



Rc = R + tq cos{ujt — 9) , 
9c = + arctan os' { t 



R + ra cos{ujt — 6) 



(70) 



We can solve Eq. (|68)) to get the coordination of the cir- 
cle's center (i?, 9), and substitute them into the position 
equation above, then we get a trajectory of drifting cir- 
cles in polar coordination. Here we study the difFusion 
of many particles, it is the circle's center that determines 
the position of particle in larger scale, so we can put the 
local circular motion aside. 

We consider a cylindrical density concentration, i.e., 
dgNs = 0, the density gradient only exist along the ra- 
dius. Further more, we withdraw the external driven 
force, Fg = 0, Fji = 0. The velocity under nonzero den- 
sity gradient is. 



R0^D, 



9R- 



R^D 



RR- 



OrNs 

Ns 



(71) 



Now we see, besides the particle flow Vr along the radius, 
there is a transverse flow Vg perpendicular to the radius. 
The transverse velocity is induced by the density gradient 
along the radius(Fig. [16] (a)). 

For a more complex density gradient distribution, the 
particle number conservation equation still holds. 



dNs 



d 



dt dR 



[NsVr) 



Rdg 



{NsVg) = 0. 



(72) 



Substituting the velocity equation ([68]) into the conser- 
vation equation above, we get the diffusion equation in 
polar coordination. 



dNs 
dt 



did d 

-^[aggNsFg + agRNsFR] 



Id, d 
RW^'"dR^' 

^^g[<^RRNsFn- 

0. 



1 d 

aRgNsFg] 



(73) 



The drift tensor and diffusion tensor are state dependent. 



Aj = A,(i?,( 



(74) 

Eq. ()73|) is an equivalent formulation of the difFusion 
equation p4p in polar coordination. 

We first consider a special case that the drift tensor 
and difFusion tensor are both constant. Then for a simple 
external driven force, 



Fg — 0, Fr = CR, 
the diffusion equation ([73]) becomes. 



(75) 



dNs 
dt 



= D 



d^Ns 



RR- 



Dg 



1 d^N, 



[D 



Re 



dR^ R^ (9612 

1 527V, 1 ON 

DgR] — — Djig 



ctrrC 



dNs 

de 



' R dOdR 
- i^brCR 



dNs 
dR 



i?2 de 

agnCNs. (76) 



We first computed the special case that C=0, this means 
there is no external force. The friction parameter is 
1] — 1. The normalization factor mA^ — 1. Angular 
frequency is uj — 10. The radial difFusion constant is 
Dji — 1. angular difFusion coefficient Dg — 0.9. The nu- 
merical value of difFusion tensors are Dru = l,Dgg = 0.9, 
DRg + DgR = —1, DRg = 9. The drift tensor are 
cee = o-flfl = 1, o-fle = -^br = 10. The particles con- 
centrate on a small half circle with radius i? = 0.1 from 
to TT. The initial value of density on the small circle is 
1 . The numerical evolution of the density distribution at 
time i = 9 is shown in Fig. [17] the density center shifted 
to TT. In an external Force field with rotational symmetry 
around origin, the trajectory oF circle particle during dif- 
fusion is not simply along the radius direction, but has 
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transverse shift perpendicular to the radius. In the re- 
gion covered by a half disc, the particle concentrating on 
the radius oi 6 = tt is much more than that on the radius 
of 61 = 0. 

We computed the velocity field distribution, one can 
see clearly the transverse flow perpendicular to ra- 
dius(Fig. [T9](a)). One special character of circle particle 
flow is there exist a transverse flow induced by velocity 
gradient. In the polar coordination, this transverse flow 
field is given by 

Je = -OnVe ^ dR[DgR— ^Dee — ——) 

- dR{aRRFR + aReFe). (77) 

The derivative of velocity is acceleration. The flow field 
J is actually the off-diagonal acceleration. The corre- 
sponding fiow field J with respect to the velocity field 
distribution Fig. [19] (a) is showed in Fig. [19] (a). The 
effective force flow field J has large projection perpendic- 
ular to radius. 
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FIG. 17: The density distribution in R~6 plane at time t = 9. 
The slop of the external force is C = 0, i.e., no external force. 

We now propose an experiment to illustrate the trans- 
verse flow of diffusing circle particles in a disk region. 
First, we collect millions of circle particles, and confine 
them on a very small disk covered the origin. Second, 
one can place two blocking wall along the radius, one is 
placed along 9 — 0, the other is placed along 9 = n. 
Then, we withdraw the confining potential in the small 
disk region, and set the particles free. If there is no block- 
ing walls, the density wave is just expanding circles with 
the origin as their common center. There is a transverse 
fiow Ve perpendicular to the radius, but we can not see 
it from density distribution, because the fiow also has ro- 
tational symmetry. The blocking wall break the radial 
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(a) (b) 

FIG. 18; (a) The density distribution in i? — ^ plane at time 
t = 9. The slop of the external force is C = 1/8. (b) The 
density distribution in R — 9 plane at time t — 9. The slop of 
the external force is C = 2.. 

symmetry, the flow is cut off when it hits the wall, thus 
we will observe the accumulation of particle on one side 
of wall, and only a few of particles on the other side of the 
wall. Fig. [TH] (b) showed the accumulation phenomena 
for particle circling in clockwise direction. If the particle 
is circling in the counterclockwise direction, the particle 
will accumulate in the opposite side of the wall. If there 
exist two kinds of particles, one is circling in clockwise 
direction, the other is circling in counter clock wise direc- 
tion, the clockwise particle will accumulate in the region 
9 = (— a, 0) and 9 — {n — a,Tr), a > 0. The counter- 
clockwise particle will accumulate in the 9 = (0,/3) and 
9={TT,n + l3), P>Q. 

The strength of external field determines whether we 
can observe the transverse flow or not. In the discussions 
above, we are studying the case of C = 0, i.e., no external 
field is applied. Now we take C = 1/8, there is a small 
external field along the radius with slop C = 1/8. We 
numerically computed the density distribution by keep- 
ing the same value of other parameters as that for C = 0. 
The amplitude of the density center on = tt is shrink- 
ing(Fig. [TO] (a)). We continue to increase the slop of 
the external field to , the density distributed in the block 
zone (0,7r) becomes almost homogeneous(Fig. [H] (b)). 
For C — 2, the external force along the radius is so strong 
that the effective transverse force can be neglected. Thus 
for a strong external force along radius, it is hard to ob- 
serve an apparent bias distribution of density. But the 
transverse fiow still exist. 



VII. THREE DIMENSIONAL DIFFUSION 
EQUATION OF MANY CIRCLING PARTICLES 

Microorganisms usually swim in three dimensional 
space, they follow a helical trajectory. Their angular ve- 
locity fluctuates according to external stimulus. A clas- 
sical electron or ion moving in three dimensions follows 
a helical trajectory. The center line of the helical path 
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of the particle can be expressed as 



(a) 



(b) 



FIG. 19: (a) The velocity field distribution in _R — 6^ plane at 
time t — 9. The slop of the external force is C = 0.(b) The 
corresponding transverse force flow field distribution in i? — 6* 
plane at time t = 9. The slop of the external force is C = 0. 

is magnetic field line. We will study a box of circling 
particles in three diniensions(Fig. [^U)) . The particle is 
swimming in thick liquid with low Reynold number. We 
focus on the most general case, and do not distinguish the 
self-propelled particles or passive swimming particles. 
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FIG. 20: (a) The circling particles swimming in three dimen- 
sions, (b) The decomposition of the motion of a circling par- 
ticle in three dimensions. The local circling motion is around 
bj. The motion of the circle's center is along the dark thick 
curve, cj is parallel to the tangent vector of the center line(the 
dark thick curve). 

We take the helical path as the combination of a local 
circling motion and a motion of the circle's center. The 
local circling motion is confined in a sphere with radius 
To . The direction of the circling motion is characterized 
by a angular velocity The circles center can move 

around in the whole space. They are steered by external 
driven force of external stimulus. In two dimensions, the 
angular velocity has only one component. The particle 
is constantly rotating in X — y plane. In three dimen- 
sions, the particle may rotates in different directions, the 
angular velocity can be projected Xo X — Y — Z three 
directions. 



(78) 



If the center of the circle is placed at the origin, and the 
angular velocity vector is Lu{t), the instantaneous position 



ro 



^z^x ~ ^x^z cos(a;i) — UyU sin(wt)]. 



ry(t) = — [wyWj; — Wa^i^j, cos(a;i) + cJxCj sin(w<)], 
rx{t) = ^[wjjCjj; + (w^ +w^)cos(a;i)]. 



(79) 



In two dimension, the particle is rotating around Z— axis, 
so the other two component of the angular velocity is 
zero, uJx — ujy — Q. Then we have uj = ui^, the posi- 
tion equation d79l) reduced exactly to the circle in two 
dimensions 

rx{t) = ro{t) cos ujt, ry{t) — ro{t) sin ujt. (80) 

We place the center of the circle to an arbitrary curve 
R{t), this curve is the trajectory of the center of the circle 
. We assume the angular velocity vector uj{t) is parallel 
to the tangent vector of the R{t). Then position of the 
circling particle is the solution of equation, 

x(t) = Rx+rx(t), 
y{t) = Ry + r-yit), 
z{t) - Rz+r,{t). (81) 

Following the same process as we model the diffusion 
phenomena in two dimensions, we took the unit step size 
to measure diffusion larger than the radius of the circle. 
Then we consider the equation of motion for the center 
line by including the diffusion force Fdi / fusion from den- 
sity gradient, friction force FfricUom external driven force 
Fexternai and Centripetal force Fcentripetai from non-zero 
angular velocity. Newton's law give us the same equation 
of motion, 



m^r— — —mVii X cj + — f— mr]VR, (82) 



dt 



except that here every vector has three components. This 
equation is composed of three coupled equations. 



mRx = 


-m{RyUJz 


- R^UJy) 


-Dx 


mRy = 


-m{RzUJx 


- Rxi^z) 


-Dy 


mRz = 


-m{RxLjy 


— RyiOx) 





dRxN, 

OrzNs 



mriRx + Fx, 
mrjRy -f- Fy, 
mriRz + F^. 



Again we assume the swim over damped, then we set the 
inertial terms as zero, 

miix = 0, mRy = 0, mR^ = 0. (83) 
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The solution of the center velocity is 



y - ^yx ^yy ^yz , 

~(" ^yx-^x ~t~ (^yyFy ^yz-^Zi 

Jl'?: ' ^ ,^ J-^ZU ^ZZ~ 



H~ ^zx-^x H" ^zyFy -\- (Tzz-^z- 



(84) 



This equation governs the motion of the center. The 
complete equation for a helical path is Eq. (|5T|) . The 
explicit formulation for the ultimate equation of a helix 
can be derived by substituting the center equation ([5^ 
and the circling equation ((7^ into Eq. (|8ip . When the 
circling velocity r ^ R, the solution of Eq. (ISTI) is a 
helical curve. 

Both the diffusion tensor and drift tensor of the center 
Eq. (|84l) are 3x3 matrix, the elements of the matrix are 



1 



[if + mul), 



Dxy — DyGxy, <Jxy — "^(^ f^f^z +^x^y), 
Dxx = DzfTxz, (Txz = -^iv^y + i^x^z), 



Dyx — Dx<Jyx, O'yx " '7^ z +^x^y), 



Dyy = DyCTyy, ^Tyy ^ + i^^) , 



Dyz = Dz<Tyz, CTyz = 7^(-?7^i + WyW^), 



1 

1 



^zx — ^x^zx, '-' zx — 7^(~'/'^l/ '^^x^z). 



CI 



Dzy ~ DyGzy, CTzy = -^{tI^x + i^y^z)i 



^ ZZ ^Z^ZZ-i ^ ZZ ^9 



and the normalization factor is 



(85) 



Cq = mr](jf' + mw^ + ?7ia;,y + muj'^). (86) 

The drifting tensor a is the sum of a symmetric matrix 
cr^ and antisymmetric matrix cr^, 



^ I rf' + mujl LUxi^y LOxi^z 



(87) 



-rjUJz T]ujy 
rjUJz -rjUx 
~r]ujy rjujx 



Following the same calculation for two dimensional dif- 
fusion theory, we can get the diffusion equation in three 
dimensions by substituting the velocity in to the mass 
conservation equation 



dtNs + Vr ■ NsV = 0. 



(88) 



If the friction coefficient 77 — ri{Rx, Ry, Rz), drift ten- 
sor aij = a{Rx,Ry,Rz) and diffusion tensor Di — 
Di{Rx, Ry, Rz) are all state-dependent, we will get a 
complicate three dimensional diffusion equation. 



d 

dR, 



0. 



(89) 



For the special case that, the friction coefficient 77, drift 
tensor aij and diffusion tensor Di are constants, we get 
an explicit expansion of the diffusion equation. 



dt 



- [Dxy + Dyx]dR^dR^Ns - Dxxdl^Ns 

- [Dxz + Dzx]dR^dRMs - Dyydl^Ns 

- [Dzy + Dyz]dR^dR^N, - Dzzdl^Ns 



[<yxxFx + O-xyFy + (JxzFz] Or^Ns 
[<7yxFx + O'yyFy + CTyzFz] dR^Ng 
[(^zxFx + CTzyFy + (TzzFz] Or^Ns 

[o-xxdR^Fx + (JxydR^Fy + axzdR^Fz] Ns 

Qr^Fx + CTyydR^Fy + CTyzdR^Fz] Ng 
[(^ZxdR^Fx + OzydR^Fy + GzzdR^Fz] Ns 

0. 



(90) 



One can choose different coordination system to refor- 
mulate the three dimensional diffusion equation for the 
convenience of physics, especially when the physical sys- 
tem has rotational symmetry. If the initial distribution 
of particle has spherical symmetry, then we can assume 
the density gradient also has spherical symmetry. Then 
we can express Eq. ((89|l by polar coordination. The den- 
sity gradient distribution looks like a spherical balloon 
expanding away from its center, but there are transverse 
flows distributed on the surfaces of the balloon. There 
is a topological constrain on the distribution of the flows 
on the surface of a sphere, it states there must exist two 
singular points at which the flows contradict one another. 
So the flow on the spherical surface like to fuse into the 
third flow along radius. 



VIII. THEORETICAL UNDERSTANDING OF 
THE SPERM DIFFUSION EXPERIMENT 

In the experiment of sperm diffusion[l[, the sperms are 
kept in a sperm reservoir in the beginning. The sperm 
reservoir is situated at one end of a rectangular migration 
channel. On the other end of the rectangular migration 
channel is a chemoattractant reservoir. By observing the 
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diffusion phenomena, they studied how the sperm be- 
haves under the action of chemoattractantm . Some other 
earher experiments showed that the chemoattractant acts 
as some effective attractive action to sperms by chemical 
stimulus. This chemoattractive force change the veloc- 
ity of sperm[l3|[l3l and makes them move to eggs. The 
experimental data of Ref. [l| shows that sperm liken to 
concentrate on one side of the channel(FiglT](a)) and the 
sperm density concentration illustrated a plateau along 
the length of channel, the plateau is near the reservoir 
source of sperm(around X2 in Fig. [T]) [l[ . 

The sperm used in experiment |l| is Arbacia punctulata 
spermatozoa. The size of the sperm head is about 5 /im, 
but its flagellum has a length of 30-60 fim. In three di- 
mensions, when there is no external stimulus, the trajec- 
tory of Arbacia punctulata spermatozoa is a helical curve. 
The center line of the helical curve is almost straight [l9| . 
When chemoattractant spreads into the solution, the 
sperm will receive stimulus, and reorient its swimming 
direction. It is believed that the sperm will bent their cir- 
cling direction to the source of chemoattractant [Tsj [l9| . 
When a sperm is swimming near a boundary, it swims in 
a circle. 
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velocity for one sperm is suppressed by the presence of 
other sperms. The sperm experiment pj study the dif- 
fusion of millions of sperms instead of one. The sperms 
in the top surface will swim in circles. The thickness 
of the top layer is about 5 fim, which is the size of one 
sperm head. The top layer behaves like a boundary for 
the sperms below, so the sperms in the second layer will 
also swim in circles. We can theoretically get at most 20 
layers of circling sperms. In reality, it is impossible pile 
up alive sperms like bricks. The sperms are not restricted 
in the top surface, they are fluctuating up and down all 
the time. So a rough estimate of the effective thickness 
of the one layer is 10 fim. If the density is high enough, 
we can have 10 layers. 

When the density is only high enough for sperms to 
form three layers, the two dimensional diffusion theory is 
still a good approximation. Suppose there are 3 layers on 
top of each other within 100 fim. The effective thickness 
of one layer is 10 /im, then we have 70 fj,m free distance. 
Suppose the distance between the nearest neighboring 
layers are equal, the distance is just 70/3 ~ 23 fim. Re- 
call that minimal square of one sperm circle is 60 x 60 
/im^, the sperms do not enough space to switch its an- 
gular velocity from ojz to {oJx,ojy). In most cases, the the 
majority of the sperms would swim in circle. The plane 
of the circling is parallel to the top surface. The two di- 
mensional diffusion theory captured the most significant 
dynamics. The fluctuation in the third dimension is a 
small perturbation. If the density is very low, we should 
use many body theory and take into account of the inter- 
action between sperms. The diffusion theory is not such 
an accurate theory for few-body system. 



FIG. 21: A cross section of the migration channel of 
experiment 0]. The length along X-axes is 7000-9000 fim. 
The width in Y-axes is 500 /xm, and the height in Z-axes is 
100 nm. The radius of the sperm circle is 30 fim. 

The rectangular migration channel in experiment [l| 
has a length of 7000-9000 //m, its width is 500 /im, and 
the height of the channel is 100 /im. This channel is a 
three dimensional channel for one sperm, but for many 
sperms, it can be well approximated by a two dimen- 
sional channel. The radius of the sperm circle is fluctu- 
ating around 30 /im. So one sperm need at least a small 
square of 60 x 60 fim^ to swim[7j 6] 17][1§|. The rectan- 
gular channel of experiment [l| has two surfaces: the top 
surface and the bottom surface. The distance between 
these two surfaces is 100 iim. A sperm head is 5 /xm, so 
there are 90 /im free distance left for the sperm to swim 
between the two surfaces. When the sperm is right in 
the middle, the distance between the sperm and one of 
the two surfaces reaches maximal, 45 fim. This distance 
is larger than the radius of the sperm circle. So the an- 
gular velocity of one sperm can rotate from Uz to uJx or 
to ijy. But when we consider two sperms, one is on top 
of another. If the two sperms try to rotate their angular 
velocity without any conflicting, the height of the chan- 
nel must be at least 120 fim. The fluctuation of angular 




FIG. 22: The density distribution Ns{Ra:,Ry,t = 15) of dif- 
fusion Eq. (|92[) for C — 2. The light color indicates high 
density, dark color indicates low density.. 

If one can reduce the height of the rectangular chan- 
nel in experiment [H from 100 fim to 5 ~ 10 fxm, then 
the angular velocity in X and Y direction, uj^jOJy, will 
be greatly suppressed. The dominant component of an- 
gular velocity is uj^. In that case, the two dimensional 
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diffusion tlieory would be a good tlieoretical approxima- 
tion. Suppose one can do tfie same experiment |1| but 
using a very tliin cfiannel with freight 5 ^ 10 /Ltm, then 
two dimensionai diffusion theory will be a good approxi- 
mation. Since there is a chemoreservior attached on the 
other end of the channel, we assume the chemo-attractive 
force distribution is 



Fx — C Rx, 



Fy=0. 



(91) 



This means the shorter distance a sperm stays away from 
chemo-attractant, the stronger chemo-attractive force it 
feels. Biologist found that almost all sperms are circling 
with same handedness [6]. So we choose angular velocity 
ujz > 0. The diffusion equation for this external driven 
force is 



dt 



DxxdlxNs + DyydlyN, 

{Dyx + Dxy)dlixdRyNs 



^ C (Jyx Rx dflyNs 

— C (JxxN_^. 



C CTxx Rx dRxN, 



(92) 



We take the diagonal diffusion coefficient as Dxx = 
0.9, Dyy = 1. The off-diagonal diffusion coefficient is 
Dyx -f Dxy = —1. The off-diagonal drifting tensor is 
<yyx — 5. The diagonal drifting tensor is Uxx — 2. The 
slop of the external force is controlled by C . From Eq. [92l 
one can see if we change C, there are three terms that will 
change simultaneously. We first study the case of C = 2. 
The numerical evolution of density distribution at time 
t — 15 is shown in Fig. [22]) . The density center shift 
to positive Y-direction. We integrate the density along 
Y-axis Na{x,t) = Ngdy, and plot the density evolu- 
tion Ns{x,t) along X-direction(Fig. [23l). The center of 
the density distribution is around Rx — 2. The diffusion 
curve along Rx is not a simple exponential decay. The 
density center can propagate as a pump. 




FIG. 23: . The density distribution of 7Vs(i?^,t = 15) = 
/J" NsdRy for diffusion Eq. ((92)) at C = 2. 

In fact, the experiment [l| measured the density dis- 
tribution for one snapshot of diffusion image. The ex- 
periment date of Ref. [l| is collected from the diffusion 
image along the length of the migration channel. A point 
on the curve represents the average density of sperms in 



each small rectangular box along Y-axis(the gray box in 
Fig. [U (a)). The data of Ref. yy] is fitted by one di- 
mensional diffusion equation. The data shows apparent 
deviation from exponential decay around 2; = 2 in Fig. [T] 
(a). They measured the diffusion of Arbacia punctulata 
spermatozoa in seawater, jelly coat solution and resact 
solution for different solutions. In all of these measure- 
ments, the high density dots deviating from exponential 
decay appears around a; = 2 in Fig. [1] (a) . Thus these 
high density dots are likely not erroneously fluctuating 
from experiment. 

The vortex generated during diffusion will reduce the 
speed of diffusion in X-direction, as particles like to flow 
around vortex core. The vortex flow results in the trap- 
ping of a large amount of particles in a small region. The 
bump in density distribution is a signal of vortex config- 
uration. The velocity of particle is given by 

Rx F)xy ~ Dxx ^~ O'xxC Rxi 



Rv - -Dyx^^r-^ ~ Dyy^^—^ + ayxC Rx. (93) 



We computed the velocity field distribution at time t = 
15 in Rx — Ry plane (Fig. [24|) . one can see a vortex 
line is generated in the high density region particles are 
flowing out of these vortex lines. These vortex lines are 
distributed along the upper boundary. The most con- 
densed particle flow is around the vortex line, it flows to 
positive i?Y-direction. 

For the external force configuration here Fx = 
C Rx, Fy = 0., we computed the transverse fiow field 
J induced by velocity gradient, 

Jy DxyUjdRxdRylnNs + Dxxi^dRxdRxlnNs 
+ ujaxxC 

Jx = Dyx^dRydRxlnNs + DyyUjdRydRylnNs. (94) 

The strong flow of J is distributed in the high density 
region, it flows to positive Ry direction(Fig. [25]) . On the 
upper boundary, all flows point to positive Rx direction. 
In the low density region, turbulent flow began come into 
being, one can see clearly a string of vortex. The flows 
around vortex snake into the flow near the upper bound- 
ary. Generally we now can understand why particles on 
the upper boundary swim faster. Experiment observed 
fast sperms on the upper boundary, they swim faster than 
others. From our numerical evolution, we see even the 
particle are all have the same angular frequency, they will 
generate a fast particle flow near the boundary. There- 
fore, it is most likely that the fast sperm is not really fast, 
they are just being pushed by other sperms, this make 
them stay in front of others. 

The data of Ref. [l| shows the highest density peak 
is right at the origin Rx = 0, which is also the position 
of the sperm reservoir. We reduced the slop of external 
fleld from C = 2 to C = 1/2, the numerical evolution of 
the diffusion equation showed a peak at the origin(Fig. 
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FIG. 24: . The velocity field distribution V = (K, Vy) corre- 
sponding to the density distribution Fig. 1221 The slop of the 
external force is C = 2. 
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FIG. 25: . The effective force flow field J — {Jx, Jy) at time 
t — 15 corresponding to the density distribution Fig. 1221 The 
slop of the external force is C = 2. 

(b)). The height of the bump is lowered down, and 
becomes a plateau (Fig. [26] (b)). Comparing Fig. [26] 
(b) with the density distribution for C — 2 (Fig. 1^51) . we 
found the tail of the bump for C — 1/2 is broadened. The 
highest density region focused on a small region in upper 
half plane(Fig. [26] (a)). Compare it with the density 
distribution for C = 2, one may see, for the case of C = 
2, the high density region covered a larger region along 
Y-direction. This means the increase of external field 
along X-axes enhanced the trapping of particles in Y- 
direction. In other words, if we increase the external 
force in X-direction, the particles do not move much in X- 
direction, instead they flow to Y-direction. This diffusion 
phenomena is special for circle swimmers. 

The data in Ref. [l| is for one snapshot of the whole 
diffusion movie. We do not have enough information to 
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FIG. 26: (a)The density distribution Ns{Rx,Ry,t = 15) of 
diffusion Eq. ([92]) for C = 1/2. The light color indicates high 
density, dark color indicates low density, (b) The density 
curve of Ns{Rx,t = 15) = /J" N^dRy for C = 1/2. 

compare theory with experiment. As far as we can see 
from the measurement of the one dimensional density 
distribution, the the highest density peak at the origin 
remains in different measurement. The sperm reservoir 
is right at the origin. Biological observation suggest that 
sperms like to stay near a boundary. I guess the ap- 
pearance of this high density peak is due to experiment 
devices. The sperms are attached by the wall of the chan- 
nel in the beginning, they did not move too much in the 
whole period of diffusion. So this peak at the origin do 
not have dynamic effect. 

Generally speaking, the density evolution of the circle 
particles diffusing across a rectangular channel is not a 
simple exponential decay, many nonlinear flow fleld are 
involved during diffusion. For a more detail study, we 
may take into account of the shape of the circling particle, 
the sperm head has a shape of rod which is beating asym- 
metrically within a small angle around the tangent vector 
of circular trajectory Q. The beating head of a boundary 
swimmer may hit the wall before its beating amplitude 
reach maximal point. As a swimmer liken to slide along 
the lengthy direction of the head, the cut off of beating 
amplitude in Y-direction by the wall will enhanced its 
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motion in i?x— direction. If a rod is surrounded by many 
rods, there exist frequent collision between rods. We can 
use Boltzmann equation to study the effect of collisions 
in diffusion process. One approximation of the linearized 
Boltzmann equation is Fokker-Planck equation[20]. The- 
oretical simulation showed that a number of sperms can 
swim together under hydrodynamic interactions between 
beating flagella 21]. So fluid dynamic theory may help us 
to understand the dynamics of many circle particles. 



IX. CONCLUSION 

The diffusion phenomena of circle particle is different 
from non-circling particles. For non-circling particle, the 
projected trajectory to X-axes is independent of its pro- 
jection to y-axes. But for circle particle, the motion 
in X-direction is coupled to the motion in Y-direction. 
There exist strong correlation between the different dif- 
fusion directions of circle swimmer. This correlation exist 
not only in two dimensions, but also in three dimensions. 
Circling motion can be theoretically modeled by non-zero 
angular velocity. We just start from a helical trajectory, 
and expressed Newton's law by local coordination of the 
circle and a global coordination of the circle's center. We 
assume the motion of the circle's center is over damped, 
so one can drop the inertial term in the equation of mo- 
tion for the circle's center. The diffusion equation of 
many circle particle is expanded by nontrivial diffusion 
tensor and drift tensor. One can actually reformulated 
this diffusion equation into Fokker-Planck equation. For 
a simple external driven force, the numerical evolution of 
the tensor diffusion suggests that the diffusion flow along 
density gradient in one direction can induce a transverse 



flow in its perpendicular direction. Recall that we did not 
distinguish weather the particle is self-propelled particle 
or not, this nonequilibrium phenomena holds both for 
lifeless particles and self-propelled particles. In a rectan- 
gular diffusion channel, the transverse flow will increase 
the density on the boundary, which in turn drive the 
swimmers near the boundary move faster. The vortex 
are generated during diffusion, we can find the vortex 
configuration both from the velocity vector field and the 
transverse flow vector field induced by velocity gradient. 
Thus the diffusion of circle particle is a highly nonlinear 
process. Comparing with the experiment measurement [l| 
of the density distribution projected to X-direction, we 
can understand the bump in the density curve [l| as a spe- 
cial phenomena for circle particles. If all circle swimmers 
have the same circling handedness, i.e., either clockwise 
or counterclockwise, when they are diffusing in a chan- 
nel along X-direction, the density center either shift to 
positive Y-direction or shift to negative Y-direction, the 
density distribution along Y-direction is not symmetric. 
The sperm experiment indeed observed sperms liken to 
concentrate to one side of the channel, and sperms swim 
faster on the boundary. We analyzed the flow field dis- 
tribution, found that the particle flow point to bound- 
ary. This can explain why sperms on the boundary swim 
faster, because they are pushed by other sperms. Fur- 
ther more, we predict one phenomena: when we mix cir- 
cle swimmers with opposite circling handedness together, 
and let them diffuse across a rectangular channel in X- 
direction, they will split to two groups, one goes to pos- 
itive Y-direction, the other goes to negative Y-direction. 
The swimmers in each group all have the same circling 
handedness. The classical circle swimmers can give us a 
vivid demonstration of plasma physics [22 1. 
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